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Abstract. We present a new algorithm for vibrational control in deuterium molecules 
that is feasible with current experimental technology. A pump mechanism is used 
to create a coherent superposition of the vibrations. A short, intense infrared 
control pulse is applied after a chosen delay time to create selective interferences. A 
'chessboard' pattern of states can be realized in which a set of even- or odd-numbered 
vibrational states can be selectively annihilated or enhanced. A technique is proposed 
for experimental realization and observation of this effect using 5 fs pulses of A = 790 
nm radiation, with intermediate intensity (5 x 10 13 W cm -2 ). 



1. Introduction 

The ability to effect quantum control at the molecular level has many potential 
applications, such as control of reaction dynamics [I], selective electron localization 
[2J, and applications in quantum information processing [3J. In order to characterize 
control mechanisms and afford useful comparison between experiment and theory the 
hydrogen diatomic (and its ion, H^) is often used as a test-bed. Such fundamental 
systems are extremely useful for identifying basic physical mechanisms because of their 
simple electronic structure. The heavier isotopes of hydrogen, such as deuterium and 
deuterium hydride are equally accessible for quantum simulations and execute slower 
nuclear motion, thus deuterium (D2) has become the preferred target for time-resolved 
studies. 

In recent years, ultrashort laser pulses have been used to initiate, image jU [5] and 
control [6j [7] the ultrafast dissociation of , with the localization of the remaining 
bound electron also manipulated [2J. However, such control of fundamental molecular 
motion has not been limited to dissociation dynamics and progress has been made on 
studies of bound wavepacket motion, with control mediated by interactions occurring 
on timescales shorter than the motion itself. For example, coherent rotational motion 
in molecules can be induced and controlled [HI E] in diatomic systems using intense-field 
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pulses on sub-rotational timescales, with the mapping of rapidly rotating (hydrogenic) 
diatomics being recently achieved [TOl [TT] . 

Several proposals for coherent control of vibration in the deuterium molecular 
ion have been put forward, with similar underlying principles [3j [121 02] • Firstly an 
ultrashort 'pump' pulse can be used to induce strong field tunnel ionization of a D 2 
target on a sub- vibrational timescale (< 25 fs), creating a D^" nuclear wavepacket in a 
coherent superposition of vibrational states. Taking t = as the moment of creation 
of the molecular ion state, the subsequent vibration can be given by the eigenfunction 
expansion: 



where, {\n)} denotes the set of vibrational eigenstates (discrete and continuous) of the 
molecular ion, and {E n } the corresponding energies. Each bound eigenstate component 
has an associated amplitude {a n } and evolves with frequency u n = E n /h~, where the 
population of each n state is \a n \ 2 . 

The evolution of the probability density is dictated by the beat frequencies {uj n — u> n >) 
between the eigenstates [H] with periods on the order of 20 - 30 fs for first order beats 
(n—n' = ±1). The anharmonicity of the potential means that the eigenstate components 
de-phase within a few vibration cycles [14] . The wavepacket will become spatially de- 
localised across the potential as time evolves but will re-form (to execute well defined 
oscillations) whenever the phases of the components in ([I]) match. That is, when the 
beats have re-phased. This quantum 'revival' (at t ~ 550 fs) was recently measured 
[T5l [IB] and well reproduced by simulations [16j [IT] . In these experiments rotational 
effects in the molecular ion were found to be negligible [16] due to the large range of 
vibrational states occupied, each with different rotational constants |18j . 

Suppose that, after some time r in the wavepacket evolution, a secondary (control) 
pulse is used to alter the vibrational wavepacket. It has been shown that an 
infrared control pulse can produce an AC Stark-effect that distorts the potential well 
(adiabatically) and hence modifies the wavepacket motion [12]. Depending on the 
wavepacket motion and position when the pulse is applied, a relative enhancement 
of lower vibrational populations may occur, i.e. vibrational cooling. While the 
adiabatic picture provides a useful physical interpretation, it does not fully explain 
the redistribution process. 

To describe this, a Raman- type scheme has been proposed j3j [121 CO] where portions 
of the wavepacket may be transferred between the lso~ g (bound) and 2po~ u (dissociative) 
potential, as seen in Figure [I] In recent simulations, high-intensity pulses (/ > 10 14 W 
cm -2 ) have been used to enforce vibrational squeezing or 'quenching' into a specific state 
[13J or to create a coherent 2 (or 3) state wavepacket [3]. These studies were conducted 
with short delay times (r < 100 fs), where the first-order vibrational beat components 
are de-phasing. 

In this article, we show that intermediate intensity (/ ~ 5 x 10 13 W cm -2 ), few-cycle 
(5 fs) control pulses have important applications for coherent manipulation. We find a 
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Figure 1. Schematic of vibrational control of (not to scale). An intense infrared 
pump laser interacts with the D2 target in its ground vibrational state. The ionization 
process creates a coherent superposition of D j" lsa g vibrational modes displaced from 
the minimum of the well. The wavepacket then oscillates, due to the displacement, 
and disperses, due to the anharmonicity of the state (silver arrow). The classical 
vibrational period (that is, the centre of mass of the oscillation) is approximately 25 
fs, while the vibrational modes will fully de-phase and re-phase ('revive') after around 
550 fs. The application of a secondary, control pulse initiates Raman transitions via 
the 2pa u surface creating Stokes (red arrow) and anti-Stokes (blue arrow) transitions. 

fascinating and hitherto unseen 'chessboard' pattern in the redistribution. In a different 
outcome to previous quenching/cooling studies, it is possible to selectively enhance 
almost exclusively even or exclusively odd numbered vibrational states. Thus the 
process relies on coherence in the energy domain, rather than wavepacket manipulation 
in the time-domain. It is found that the mechanism of this process relies on 
constructive/destructive interference. This effect can be exploited at a fractional revival 
time to create strong-contrast interference patterns. 

2. Simulations of Vibrational Control 

In this article, we simulate the control of a coherent vibrational wavepacket using 
a 5 fs, 5xl0 13 W cm -2 pulse. Referring to Figure 1, an ultrashort 'pump pulse' (in 
this case provided by a Ti: Sapphire laser) interacts with the initial D 2 target in its 
ground vibrational state. The ionization process creates a coherent superposition of 
lscx ff vibrational states, as given in ([I]), displaced from the minimum of the potential 
well. The wavepacket then oscillates, due to the displacement, and disperses due to 
the anharmonicity of the state (silver arrow). The application of a secondary, 
control pulse initiates Raman transitions via the 2pa u surface (blue/red arrows), acting 
to redistribute the vibrational distribution of the wavepacket. This new vibrational 
distribution is studied as a function of the control pulse delay time, r. 

The initial D 2 target can be assumed to reside in the ground vibrational state of 
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the X l Yt g potential. In the Franck-Condon (FC) approximation, the ionization event 
would simply project this vibration function onto the manifold of vibrational states of 
the ion. Regarding this pump process, it is well known (see [191 ED]) f° r example) 
that the tunnelling-ionization rate varies with bond-length, even for the compact 
wavepacket presented by the D 2 ground vibrational state. This can modify the simple FC 
approximation that assumes the wavepacket is projected without modulation. We have 
carried out simulations using modulated wavepackets (displaced/broadened/skewed) as 
well as the FC wavepacket and find that the 'chessboard' interference effect we report 
here can be produced for each distribution. The FC approximation will be used here to 
demonstrate the 'chessboard' effect but it is not a necessary condition. A wavepacket 
with a range of populated vibrational states (such as those produced by intense infrared 
lasers) is sufficient. For initial distributions in a narrow range of vibrational states, the 
interference effects may not extend over as many states but will still occur. Indeed, for 
any experimental investigation the pump pulse parameters should be carefully chosen 
and compared to simulations for appropriate initial wavepacket conditions. 

After the pump process, the FC wavepacket resides in the ground electronic state 
of the ion. At intermediate intensities, it is sufficient to solve the time-dependent 
equation for evolution on the two lowest potential curves (lsa g and 2pa u ) within 
the Born-Oppenheimer approximation. With the infrared frequencies considered here, 
these two low lying levels need only be considered as higher-lying electronic states are 
inaccessible. Thus the configuration space can be partitioned into g and u components. 
Then the Hamiltonian, H, can be partitioned into vibration, electronic motion, and 
laser interaction as follows: 

H = T R + H e + V . (2) 

This has the matrix representation: 

b=\- : i + i £ ' (iS) A + ( ° md{R) I , (3) 




e u (R) J V F(t)d(R) 

where T R is the vibrational kinetic energy operator, e g (R) and e u (R) are the potential 
energy curves for the ls<r g and 2pa u electronic states respectively, and d(R) is the 
dipole moment between these states [21]. The external electric field F(t) is created by 
a Ti: Sapphire laser with A = 790 nm, modulated by a Gaussian profile centred at t = r 
(the delay time), with full-width-half-maximum of W (the duration), that is: 

F(t) = F cos(cj(t-r))exp{-21n2(t-r) 2 /iy 2 } . (4) 

The (cycle-average) intensity, /, of the pulse is related to the electric field amplitude, 
F , by the formula, I = ^€qcFq. The evolution of the equations is computed by a 
symmetric, split-step algorithm [22] : 

\ip{t + At)) = exp(-iT R At/2) exp(-iF e At/2) exp(-iV(t) At) (5) 
x exp(-iF e At/2) exp(-iT R At/2) \ip(t)) + 0(At 2 ) . 
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The error term in this expression arises partly from the splitting (factorization) and 
from assuming the Hamiltonian changes slowly over the time step (time-ordering error). 
In practice, the time-step At is chosen sufficiently small so that, At(dH / dt) <C H . 

The splitting and the use of a uniform grid means the highly-efficient fast Fourier 
transform [23] can be employed. The electronic coupling term V is diagonal in the radial 
dimension, and a diagonalization of the 2x2 submatrix gives an efficient scheme for 
propagation: 

F(t)d(R) 
F{t)d{R) 




exp(-iF(t)d(R)At) 

exp(iF(t)d(R)At) 




(6) 



where 1 is the unit matrix. The vibrational populations are projections of the g-state 
wavepacket on the manifold, {|^)}, while the dissociation yield will be defined by the 
population in the u state along with the g-state continuum. 



3. The Chessboard 



The vibrational populations, following the application of a short control pulse (W = 5 
fs), with I — 5 x 10 13 W cm -2 , were calculated for a range of delay times fs < r < 650 
fs. The results are presented in Figure [2j The colour density in Figure [2] (a) represents 
the final population of a vibrational level with respect to the control pulse delay. The 
distributions resulting from r = 293 fs and r = 306 fs are extracted from Figure [2] (a) 
and displayed as the bar charts [2] (c) and [2] (d), respectively. The initial probability 
distribution is shown in Figure [2] (b). 

In previous studies, the modulations in population had been attributed to the 
classical vibrational period of each level [13] . Here a description for the mechanism will 
be given in terms of interference processes in the energy domain regulated by the beat 
frequencies. 

For short delays (r < 50 fs) the anharmonic correction is small, and dephasing is not 
significant. Under such conditions, the control pulse spreads the population evenly with 
only slight enhancements/deficits. In the region r ~ 100 fs, enhancement is concentrated 
in a specific level (e.g. n = 4 at r = 100 fs). These effects are repeated around the 
wavepacket revival (~ 550 fs), were the wavepacket begins to re-form and then de-phase 
again. Higher-intensity control pulses can be used to optimize state-selective control 
[3J H3] and drive population towards one or more specific levels. 

The most striking feature of Figure [2] (a) is the alternating light and dark squares, 
suggestive of a 'chessboard' pattern, centred around r ~ 300 fs. At this fractional 
revival time, the even-numbered vibrational states will be in phase with each other but 
in anti-phase with the 'odd' states. This phase relationship can be exploited to produce 
interference effects. The effect is illustrated in Figure [2] (c), where at r = 293 fs the 
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Figure 2. Vibrational distributions following aA = 790 nm, W = 5 fs, I = 5 X 10 13 
W cm -2 control pulse, as a function of delay time r. (a) The colour density represents 
the relative population of each level, (b) Initial probability distribution. For short 
delay times, r < 100 fs, and near the revival, r ~ 550 fs, the control pulse redistributes 
broadly with only slight enhancement in any single level. In contrast, the 'chessboard' 
pattern (centred around the fractional-revival t ~ 300 fs) shows an interference pattern 
in which even or odd numbered states are annihilated. For r = 293 fs and r = 306 fs, 
cuts through the colourmap are shown in (c) and (d). 



non-neglible populations are exclusively in even numbered states. Likewise at r = 306 
fs only odd numbered states are populated, with n — 1, 3, 5, 7, strongly favoured. 

The effect can be understood by the destructive/constructive interference between 
nearest-neighbour states. Suppose we consider the effect of the pulse on isolated states 
n = 3,4,5, separately. In Figure [3] the final amplitudes and phases for each state are 
illustrated by bar charts and clocks respectively. The clocks display the phase shifts for 
the wavepacket population moved into each level, with the relative phase shift defined 
with respect the expected phase for an unperturbed state. All clocks will remain at 12 
o'clock in the absence of any perturbation. In each case, starting with a single vibrational 
state, the final probability distribution is independent of r. We also note a selection rule 
favouring neighbouring levels. To illustrate this, let us focus on the state n = 4. For the 
reader's convenience, the bar charts showing the contributions into |4) are highlighted 
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Figure 3. Vibrational transfer starting from (a) n = 3, (b) n = 4 and (c) n = 5 by 
a 5 fs, A = 790 nm, I = 5 X 10 13 W cm" 2 control pulse. The bar charts represent the 
probabilities \a n \ 2 while the clocks display the phase of the coefficients a n after the 
pulse for t = 293 fs and r = 306 fs. The distributions are concentrated around the 
initial state, while the initial state is phase-advanced. This illustrates that a coherent 
sum of the final distributions will depend strongly on r. For the reader's convenience, 
the contributions into |4) are highlighted in red, see text for details. 



in red in Figure [3] At the delay time r = 293 fs, the contributions to |4) from all three 
initial states have roughly the same phase, i.e. they interfere constructively. In contrast, 
at r = 306 fs, the contributions from |3) and |5) are almost opposite to the contribution 
from |4), giving a destructive interference. For the final population in |5), the situation 
is opposite, yielding destructive and constructive interference for r = 293 fs and 306 fs 
respectively. This pattern is repeated so that even numbered states get enhanced and 
odd numbered states get quenched for r = 293 fs, and vice-versa for r = 306 fs giving 
rise to the chessboard effect seen in Figure [2} 

Some previous work on wavepacket engineering focussed on manipulation in 
configuration space. The chessboard is an effect that pertains to engineering in the 
energy dimension. We note that, although the chessboard state is highly regular in 
energy space, the wavepacket in configuration space is still quite irregular as time 
progresses. 



4. Model 



A simple mathematical model of the interference can be obtained from perturbation 
theory. As a preliminary comment, the hierarchy of timescales in this system are worth 
noting, from the fast electronic motion to the medium-speed optical field, through to 
the comparatively 'slow' vibrational motion. The 'electronic time' is of the order of the 
dissociation energy of the molecule, with an associated angular frequency: ~ 25 x 10 15 
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rad s" 1 ~ 0.6 a.u., while the optical cycle time is uj opt ~ 2.4 x 10 15 rad s" 1 ~ 0.06 a.u., 
and finally the vibration time u e ~ 0.3 x 10 15 rad s _1 ~ 0.007 a.u. Consequently the 
pulse interaction is sudden compared to the vibrational timescale, but adiabatic in the 
context of the electronic motion. 

The effect of a laser pulse after a delay r can now be considered. A simple 
mathematical model can be constructed where the the modulated sinusoidal signal of the 
laser pulse is considered as a sequence of alternating square waveforms. The expansion 

m)) = Y. a ^y^ nt \i^) , (?) 

7,ra 

can be proposed as the solution to 

i| i^(*)> = H( t ) mt)) , (8) 

where 7 G {g, u} is the electronic state with corresponding vibration state (bound and/or 
dissociative) n. Equation ([8l can be written as a set of integral equations (t > to): 



Q 7,rc(0 = Q 7,n(fr)) - I ^ / di / a 7 / )Tl /(t / )K r) n;7',n'(* / ) ex P[ i ^7,^7 / ,n'* / ] , (9) 
7 ',n' 

where A 7in;7 / jn / = E lyH — Eji, n ' an d the coupling potential V is expanded in the basis set 
in equation ([7]). In the two-state approximation the only non-zero element has the form 

V g>n;u , k (t') = V Ujk;g>n (t') = F(t)(k\d(R)\n) , (10) 

where \k) is an energy- normalised continuum function of the it-state, with wavenumber 
k, and \n), is a vibrational eigenstate (bound or continuum) of the g-well. 

The photon energy (~ 0.06 a.u. in this case) is below the single-photon resonant 
dissociation for all states of the g potential that are shown in Figure |2j and for laser 
intensities lower than 10 14 W cm -2 , the upper electronic level it acts as a virtual state. 
The expansion coefficients of ([7]) of the (/-state can then be calculated to second-order 
accuracy. 

/T+W 
dt' exp [iA fljn;Ujfe t'] V g , n . u<k (t') (11) 

X J dfexp [\A Ut k;g,n't"\ V Ut k;g,n'{t')a g ,n<{ti) ■ 

For simplicity, consider the pulse to be applied impulsively after a delay r, for a 
duration W (corresponding to half an optical cycle), with amplitude F . A further 
simplification can now be made in the summation over the intermediate states. The 
'closure approximation' replaces the spectrum of intermediate energies by a constant, 
E u< k = E, representative of the dominant (resonant) channels. In fact the choice of E 
can be used as a free parameter if desired, but it must lie in the continuum. In the spirit 
of the approximation, we choose E as the threshold of dissociation in the u state, also 
reflecting the enhanced density of states at threshold. Then taking this as the reference 
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Figure 4. Second-order coupling strength between vibrational cigcnstatcs during 
the laser pulse (d 2 )ij. Notice the rapid decrease in coupling strength as \i — j\ grows 
leading to an approximate selecton rule i — j 6 {—1, 0, 1} for intermediate fields This 
arises from the linear form of the function d(R). Also notice that the coupling strength 
increases with vibrational levels, giving a rise to a slight heating, i.e. propensity towards 
higher vibrational levels. 



point for the zero of energy, we simply set E = 0. The closure (completeness) relation 
can then be applied so that: 

^2(n\d(R)\k)(k\d(R)\n') = {n\d\R)\n') ee < n , . (12) 



The index g may now be dropped for convenience and we have: 
ln ( r + W') = a n (0) - MoXn'-Po / dt'e 1 ^' I dt"e- iE ^" 



(13) 



and this may be evaluated without further approximation. 

As seen in Figure 4j the matrix element d 2 nn „ provides an approximate selection 
rule: n' = n, n ± 1, so that the expression can be reduced even further. Moreover, since 
E n W 3> 1, then the solution can be expressed in the simple form: 

a n (r + W) w a n (0) (1 - m nn ) - m^a^O)^-^- 1 ^ (14) 
-m n ,„ +1 a„ +1 (0)e i(£; "-^+ l)T , 

where, 

'^o<n'\ fexp{i(E n - E n ,)W] - 1 



In the case of a very short impulse, corresponding to a half-cycle, (E n — E n ^i)W 1, 
and: 



F 2 W'\ 2 

r? I a n,n' ■ 



(16) 



The term n n>n simply reflects the change in phase due to the quadratic Stark effect. 
Since this is a second-order effect, the sign of F has no consequence, and K n n < 0, 
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and hence there is a phase advance. This is illustrated in Figure [3] where the parent 
state is phase shifted forward, with the excited state n = 5 more advanced than the 
lower-energy state n = 3, in agreement with this expression for the energy dependence. 
We note that, in this model, n ntn > is not sensitive to the sign of F , that is, to the carrier 
envelope phase. Thus a sequence of alternating half-cycles creates the same effect as a 
full-wave rectified pulse. 

Consider the effect of one half-cycle. The equations above have the simple 
explanation that a n (r + W) is created by an interference pattern with its immediate 
neighbours. This pattern can be regular if the terms on the right hand side are 
coherent. In order to create strong destructive interference, the terms (K nn _ 1 a n _ 1 (0) + 
K„ in+ ia n+ i(0)) must be approximately the same size as the remaining term a„(0)(l— 
for n G 2 ... 6. This gives an optimal value of F and W. Next, the phases of the two 
interference terms must be equal, which implies the condition (for all n): 

(E n -E n ^)T = -(E n+1 -E n )r , mod(2vr) . (17) 

A rough estimate of this condition on r, can be made by the anharmonic expansion, 
E n « — D e + Huj e (n + — fajj e x e {n + |) 2 , where D e is the dissociation energy. This is 
satisfied during the fractional revival around r « -K/{2uj e x e ) ~ 280 fs. Finally, these 
terms must be out of phase with a n to create destructive interference. The Stark effect 
(14) creates a phase shift, — K n>n , for the state n. Thus in order to get destructive 
interference, the requirement is then 

K n , n = (E n - S n _i)r + vr/2 , mod(27r) . (18) 



If destructive interference occurs for even n, then according to equation (17), 
constructive interference will be observed for odd n. Furthermore, if destructive 
interference is observed for a state at r = r , constructive interference will be observed 
for r = To + 7r/co> e « tq + 11 fs. This is a qualitative explanation, and is limited by the 
crudeness of the closure approximation. Further, the assumptions we make for the pulse 
are simple. Nevertheless, it appears that this simple model explains the main features 
that we observe in the numerical simulations. 



5. Proposed Experimental Technique for Observing the Chessboard 

Having simulated the chessboard pattern and explained the underlying mechanism, 
attention can be turned to techniques for experimental observation of this fascinating 
effect. In recent experiments HB] . coherent vibrational wavepacket motion in 
has been initiated and imaged in a pump-probe configuration using short intense-field 
pulses. In these studies the vibrational wavepacket motion, and notably the vibrational 
revival, has been imaged by using a probe pulse to enforce photodissociation or Coulomb 
explosion of the vibrating molecule across a range of probe delay times. 

Fourier analysis of the fragmentation yield can return a measure of both the 
temporal and spatial nature of the wavepacket motion, which can ultimately lead to full 
characterization of the molecular motion [24). In particular, the beat frequencies that 
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Figure 5. Spectral density for the dissociation yield correlation function. The 
coherent wavepacket is controlled by a 5 fs, A = 790 nm, I = 5 X 10 13 W cm~ 2 
pulse at r — 293 fs and r = 306 fs for plots (a) and (b) respectively. The probe pulse is, 
A = 790 nm, I = 4 X 10 14 W cm~ 2 , duration 5 fs. Fourier analysis of the dissociation 
signal returns the vibrational beats uj n +2 — <-^>n present in the wavepacket motion in 
each case. The probe delay time ranges from 310 fs to 4000 fs. A shorter range of 
times, 310 fs to 1000 fs for example, gives similar features. In this case the peaks are 
broadened but still well resolved. 



dictate the motion of the probability density can be extracted from the photodissociation 
signal. 

Thus it may be possible to conduct a pump-control-probe experiment where a new 
coherent vibrational distribution (created at a specific control pulse delay time) may be 
monitored by recording the total dissociation yield as a function of probe delay time 
(V) with subsequent spectral analysis of the signal. 

It is important that for modelling any such experiment, the pump process is well 
characterised. For pulse intensities of ~ lxlO 14 W cm -2 the independence (see [25] 
and references therein) of the pump process has been considered in recent studies of 
wavepackets [13], [Ml EE] , and indeed in some recent experiments [15] the results appear 
to deviate from an FC wavepacket. It is noted here that the R-dependence of the pump 
process is sensitive to pulse intensity [27] and, as pulse intensity increases, the resulting 
ion may tend towards an FC distribution [27], However, if the intensity is set too high 
there is a probability that direct fragmentation of the D 2 target may occur in the initial 
pump pulse, clouding the detection of any bound wavepacket dynamics from the control 
or probe pulse interactions. 

At this point, it is instructive to observe that in other recent experiments [161 E] 
with pump pulses (A = 790 nm) of duration 10 — 15 fs, and intensities of 5 — 8 x 10 14 
W cm -2 the wavepacket dynamics are in fact well reproduced by Franck-Condon 
simulations. In these experiments the problem of direct D 2 fragmentation from an 
intense pump pulse is overcome by orienting the pump pulse polarisation perpendicular 
to the detection axis, ensuring that any direct fragmentation arising from this pulse goes 
undetected. 
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In this context, the proposed experimental technique for identifying the chessboard 
effect has been simulated for an FC wavepacket. This was done by fixing the control 
pulse at a desired r value, and introducing an intense probe pulse (5 fs duration, 
I = 4 x 10 14 W cm" 2 ) at a variable delay time t'. The dissociation yield was deduced 
by subtracting the remaining bound wavepacket population from the initial norm of 
the wavepacket. The Fourier transform (energy spectral density) of the dissociation 
signal is shown in Figure [5] for control delay times of (a) r = 293 fs and (b) r = 306 
fs. The calculations represent the signal correlation over very long times (4000 fs). We 
have performed the same analysis over a much shorter time (1000 fs) and, although the 
features are broadened, the signatures are still clearly visible. 

The individual states in the vibrational distribution are not resolved in isolation but 
rather the peaks in Figure [5] (a) correspond to the beat frequencies between even states. 
These are notably absent in Figure [5] (b), where only the odd numbered eigenstate beat 
frequencies are observed, in keeping with the expected vibrational distributions (from 
Figure pi). The amplitude of the beat Lu n+2 — uj n in the wavepacket motion is given by 
|a n a n+ 2| • The Fourier transform technique will not extract the absolute population 
values, |a n a n+ 2| but rather, will provide insight into which states are significantly 
populated. This is a useful technique for experimentally identifying the preferential 
creation of either even or odd numbered states. It should be noted however that periods 
of the uj n+ 2 — uj n beats are typically 10 — 12 fs and thus for any experimental verification, 
it is imperative that pulses as short as 5 fs are used. 

In principle, this technique using the photodissociation signal could be a useful 
method for identifying the main components of the wavepacket, and with state-of-the- 
art pulses now available for durations less than 10 fs, this should be achievable in the 
laboratory in the not too distant future. Indeed, with continuing advances in intense- 
field laser-pulse experiments, such an even (or odd) state wavepacket (evolving with ~ 
12 fs oscillations) may even be characterized by using a high-intensity probe pulse and a 
very recently reported 'time-series analysis' of the Coulomb explosion signal [26], where 
full reconstruction of the wavepacket has been proposed. 

6. Conclusion 

In this article, the application of an intermediate intensity (5 x 10 13 W cm -2 ) few-cycle 
(5 fs) pulse to control a coherent vibrational wavepacket in has been investigated. 
The redistribution of vibrational levels proceeds by a Raman transfer between near- 
neighbour levels. This creates interferences such that it is possible to produce strong- 
contrast fringes of the even/odd numbered vibrational levels. The effect is a property of 
the molecular ion and does not require a highly-selective ionization process. All that is 
required is that a band of vibration states is populated. Simulations using the Franck- 
Condon approximation are shown for illustration but this is not a necessary condition. 

The 'chessboard' pattern that emerges around r ~ 275 — 300 fs provides a 
pathway for creating almost exclusively even or exclusively odd numbered vibrational 
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distributions in a coherent wavepacket. Thus the process relies on coherence 
in the energy domain, rather than wavepacket manipulation in the time-domain. 
The subsequent propagation of such a created wavepacket may be probed via 
photodissociation and here it has been shown that spectral analysis of the dissociation 
data may serve to identify the populated states. This provides a pathway for using 
currently available ultrashort laser pulse technology to investigate this 'chessboard' 
effect. 
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